sOCP: a framework predicting smORF coding potential based on TIS and in-frame features and effectively applied in the human genome.

Small open reading frames (smORFs) have been acknowledged to play various roles on essential biological pathways and affect human beings from diabetes to tumorigenesis. Predicting smORFs in silico is quite a prerequisite for processing the omics data. Here, we proposed the smORF-coding-potential-predicting framework, sOCP, which provides functions to construct a model for predicting novel smORFs in some species. The sOCP model constructed in human was based on in-frame features and the nucleotide bias around the start codon, and the small feature subset was proved to be competent enough and avoid overfitting problems for complicated models. It showed more advanced prediction metrics than previous methods and could correlate closely with experimental evidence in a heterogeneous dataset. The model was applied to Rattus norvegicus and exhibited satisfactory performance. We then scanned smORFs with ATG and non-ATG start codons from the human genome and generated a database containing about a million novel smORFs with coding potential. Around 72 000 smORFs are located on the lncRNA regions of the genome. The smORF-encoded peptides may be involved in biological pathways rare for canonical proteins, including glucocorticoid catabolic process and the prokaryotic defense system. Our work provides a model and database for human smORF investigation and a convenient tool for further smORF prediction in other species.

Molecular characterization of the complete genome of a novel ormycovirus infecting the ectomycorrhizal fungus Hortiboletus rubellus.

Advancements in high-throughput sequencing and the development of new bioinformatics tools for large-scale data analysis play a crucial role in uncovering virus diversity and enhancing our understanding of virus evolution. The discovery of the ormycovirus clades, a group of RNA viruses that are phylogenetically distinct from all known Riboviria members and are found in fungi, highlights the value of these tools for the discovery of novel viruses. The aim of this study was to examine viral populations in fungal hosts to gain insights into the diversity, evolution, and classification of these viruses. Here, we report the molecular characterization of a newly discovered ormycovirus, which we have named "Hortiboletus rubellus ormycovirus 1" (HrOMV1), that was found in the ectomycorrhizal fungus Hortiboletus rubellus. The bipartite genome of HrOMV1, whose nucleotide sequence was determined by HTS and RLM-RACE, consists of two RNA segments (RNA1 and RNA2) that exhibit similarity to those of previously studied ormycoviruses in their organization and the proteins they encode. The presence of upstream, in-frame AUG triplets in the 5' termini of both RNA segments suggests that HrOMV1, like certain other ormycoviruses, employs a non-canonical translation initiation strategy. Phylogenetic analysis showed that HrOMV1 is positioned within the gammaormycovirus clade. Its putative RNA-dependent RNA polymerase (RdRp) exhibits sequence similarity to those of other gammaormycovirus members, the most similarity to that of Termitomyces ormycovirus 1, with 33.05% sequence identity. This protein was found to contain conserved motifs that are crucial for RNA replication, including the distinctive GDQ catalytic triad observed in gammaormycovirus RdRps. The results of this study underscore the significance of investigating the ecological role of mycoviruses in mycorrhizal fungi. This is the first report of an ormycovirus infecting a member of the ectomycorrhizal genus Hortiboletus.

Full-length transcriptome analysis of a bloom-forming dinoflagellate Prorocentrum shikokuense (Dinophyceae).

Prorocentrum shikokuense (formerly P. donghaiense) is a pivotal dinoflagellate species associating with the HABs in the East China Sea. The complexity of its large nuclear genome hindered us from understanding its genomic characteristics. Full-length transcriptome sequencing offers a practical solution to decipher the physiological mechanisms of a species without the reference genome. In this study, we employed single-molecule real-time (SMRT) sequencing technology to sequence the full-length transcriptome of Prorocentrum shikokuense. We successfully generated 41.73 Gb of clean SMRT sequencing reads and isolated 105,249 non-redundant full-length non-chimeric reads. Our trial has led to the identification of 11,917 long non-coding RNA transcripts, 514 alternative splicing events, 437 putative transcription factor genes from 17 TF gene families, and 34,723 simple sequence repeats. Additionally, a total of 78,265 open reading frames were identified, of them 15,501 were the protein coding sequences. This dataset is valuable for annotating P. shikokuense genome, and will contribute significantly to the in-depth studies on the molecular mechanisms underlining the dinoflagellate bloom formation.

UPF1 regulates mRNA stability by sensing poorly translated coding sequences.

Post-transcriptional mRNA regulation shapes gene expression, yet how cis-elements and mRNA translation interface to regulate mRNA stability is poorly understood. We find that the strength of translation initiation, upstream open reading frame (uORF) content, codon optimality, AU-rich elements, microRNA binding sites, and open reading frame (ORF) length function combinatorially to regulate mRNA stability. Machine-learning analysis identifies ORF length as the most important conserved feature regulating mRNA decay. We find that Upf1 binds poorly translated and untranslated ORFs, which are associated with a higher decay rate, including mRNAs with uORFs and those with exposed ORFs after stop codons. Our study emphasizes Upf1's converging role in surveilling mRNAs with exposed ORFs that are poorly translated, such as mRNAs with long ORFs, ORF-like 3' UTRs, and mRNAs containing uORFs. We propose that Upf1 regulation of poorly/untranslated ORFs provides a unifying mechanism of surveillance in regulating mRNA stability and homeostasis in an exon-junction complex (EJC)-independent nonsense-mediated decay (NMD) pathway that we term ORF-mediated decay (OMD).

Characterization of a novel gammapartitivirus infecting the phytopathogenic fungus Pyricularia oryzae.

In this study, we identified a novel double-strand RNA (dsRNA) mycovirus in Pyricularia oryzae, designated "Magnaporthe oryzae partitivirus 4" (MoPV4). The genome of MoPV4 consists of a dsRNA-1 segment encoding an RNA-dependent RNA polymerase (RdRP) and a dsRNA-2 segment encoding a capsid protein (CP). Phylogenetic analysis indicated that MoPV4 belongs to the genus Gammapartitivirus within family Partitiviridae. The particles of MoPV4 are isometric with a diameter of about 32.4 nm. Three-dimensional structure predictions indicated that the RdRP of MoPV4 forms a classical right-handed conformation, while the CP has a reclining-V shape.

Evolution of a Human-Specific De Novo Open Reading Frame and Its Linked Transcriptional Silencer.

In the human genome, two short open reading frames (ORFs) separated by a transcriptional silencer and a small intervening sequence stem from the gene SMIM45. The two ORFs show different translational characteristics, and they also show divergent patterns of evolutionary development. The studies presented here describe the evolution of the components of SMIM45. One ORF consists of an ultra-conserved 68 amino acid (aa) sequence, whose origins can be traced beyond the evolutionary age of divergence of the elephant shark, ~462 MYA. The silencer also has ancient origins, but it has a complex and divergent pattern of evolutionary formation, as it overlaps both at the 68 aa ORF and the intervening sequence. The other ORF consists of 107 aa. It develops during primate evolution but is found to originate de novo from an ancestral non-coding genomic region with root origins within the Afrothere clade of placental mammals, whose evolutionary age of divergence is ~99 MYA. The formation of the complete 107 aa ORF during primate evolution is outlined, whereby sequence development is found to occur through biased mutations, with disruptive random mutations that also occur but lead to a dead-end. The 107 aa ORF is of particular significance, as there is evidence to suggest it is a protein that may function in human brain development. Its evolutionary formation presents a view of a human-specific ORF and its linked silencer that were predetermined in non-primate ancestral species. The genomic position of the silencer offers interesting possibilities for the regulation of transcription of the 107 aa ORF. A hypothesis is presented with respect to possible spatiotemporal expression of the 107 aa ORF in embryonic tissues.

ICTV Virus Taxonomy Profile: Hantaviridae 2024.

Hantaviridae is a family for negative-sense RNA viruses with genomes of about 10.5-14.6 kb. These viruses are maintained in and/or transmitted by fish, reptiles, and mammals. Several orthohantaviruses can infect humans, causing mild, severe, and sometimes-fatal diseases. Hantavirids produce enveloped virions containing three single-stranded RNA segments with open reading frames that encode a nucleoprotein (N), a glycoprotein precursor (GPC), and a large (L) protein containing an RNA-directed RNA polymerase (RdRP) domain. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Hantaviridae, which is available at ictv.global/report/hantaviridae.

ICTV Virus Taxonomy Profile: Fusariviridae 2024.

Fusariviridae is a family of mono-segmented, positive-sense RNA viruses with genome sizes of 5.9-10.7 kb. Most genomic RNAs are bicistronic, but exceptions have up to four predicted ORFs. In bicistronic genomes, the 5'-proximal ORF codes for a single protein with both RNA-directed RNA polymerase (RdRP) and RNA helicase (Hel) domains; little is known about the protein encoded by the second ORF. Fusarivirids do not appear to form virions. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Fusariviridae, which is available at ictv.global/report/fusariviridae.

Biological characterization and complete genome analysis of the newly isolated Serratia liquefaciens phage vB_SlqS_ZDD2.

A novel lytic phage named vB_SlqS_ZDD2 was isolated from hospital sewage using the double-layer agar method with Serratia liquefaciens ATCC 27592 as the host. BLASTn analysis showed that the genome sequence of phage vB_SlqS_ZDD2 did not resemble any other phages in the NCBI database. Phenotype and phylogeny analysis indicated that this phage might be a new member of the class Caudoviricetes. Phage vB_SlqS_ZDD2 has a dsDNA genome of 49,178 bp with 55% GC content and has 73 open reading frames. This phage exhibited strong lytic activity and a wide range of pH (3-12) and temperature tolerance (below 70℃).

Complete genome analysis of a novel narnavirus in sweet viburnum (Viburnum odoratissimum).

Trees and shrubs provide important ecological services. However, few studies have surveyed the virome in trees and shrubs. In this study, we discovered a new positive-sense RNA virus originating from Viburnum odoratissimum, which we named "Vo narna-like virus". The complete genome of Vo narna-like virus is 3,451 nt in length with an open reading frame (ORF) encoding the RNA-dependent RNA polymerase (RdRP) protein. Phylogenetic analysis placed this virus within the betanarnavirus clade, sharing 53.63% amino acid sequence identity with its closest relative, Qingdao RNA virus 2. The complete sequence of the virus was confirmed by rapid amplification of cDNA ends (RACE) and Sanger sequencing. Small interfering RNA (siRNA) analysis indicated that this virus interacts with the RNA interference (RNAi) pathway of V. odoratissimum. This is the first report of a narnavirus in V. odoratissimum.

Widespread stable noncanonical peptides identified by integrated analyses of ribosome profiling and ORF features.

Studies have revealed dozens of functional peptides in putative 'noncoding' regions and raised the question of how many proteins are encoded by noncanonical open reading frames (ORFs). Here, we comprehensively annotate genome-wide translated ORFs across five eukaryotes (human, mouse, zebrafish, worm, and yeast) by analyzing ribosome profiling data. We develop a logistic regression model named PepScore based on ORF features (expected length, encoded domain, and conservation) to calculate the probability that the encoded peptide is stable in humans. Systematic ectopic expression validates PepScore and shows that stable complex-associating microproteins can be encoded in 5'/3' untranslated regions and overlapping coding regions of mRNAs besides annotated noncoding RNAs. Stable noncanonical proteins follow conventional rules and localize to different subcellular compartments. Inhibition of proteasomal/lysosomal degradation pathways can stabilize some peptides especially those with moderate PepScores, but cannot rescue the expression of short ones with low PepScores suggesting they are directly degraded by cellular proteases. The majority of human noncanonical peptides with high PepScores show longer lengths but low conservation across species/mammals, and hundreds contain trait-associated genetic variants. Our study presents a statistical framework to identify stable noncanonical peptides in the genome and provides a valuable resource for functional characterization of noncanonical translation during development and disease.

cirCodAn: A GHMM-based tool for accurate prediction of coding regions in circRNA.

Studies focusing on characterizing circRNAs with the potential to translate into peptides are quickly advancing. It is helping to elucidate the roles played by circRNAs in several biological processes, especially in the emergence and development of diseases. While various tools are accessible for predicting coding regions within linear sequences, none have demonstrated accurate open reading frame detection in circular sequences, such as circRNAs. Here, we present cirCodAn, a novel tool designed to predict coding regions in circRNAs. We evaluated the performance of cirCodAn using datasets of circRNAs with strong translation evidence and showed that cirCodAn outperformed the other tools available to perform a similar task. Our findings demonstrate the applicability of cirCodAn to identify coding regions in circRNAs, which reveals the potential of use of cirCodAn in future research focusing on elucidating the biological roles of circRNAs and their encoded proteins. cirCodAn is freely available at https://github.com/denilsonfbar/cirCodAn.

Stem-loop-induced ribosome queuing in the uORF2/ATF4 overlap fine-tunes stress-induced human ATF4 translational control.

Activating transcription factor 4 (ATF4) is a master transcriptional regulator of the integrated stress response, leading cells toward adaptation or death. ATF4's induction under stress was thought to be due to delayed translation reinitiation, where the reinitiation-permissive upstream open reading frame 1 (uORF1) plays a key role. Accumulating evidence challenging this mechanism as the sole source of ATF4 translation control prompted us to investigate additional regulatory routes. We identified a highly conserved stem-loop in the uORF2/ATF4 overlap, immediately preceded by a near-cognate CUG, which introduces another layer of regulation in the form of ribosome queuing. These elements explain how the inhibitory uORF2 can be translated under stress, confirming prior observations but contradicting the original regulatory model. We also identified two highly conserved, potentially modified adenines performing antagonistic roles. Finally, we demonstrated that the canonical ATF4 translation start site is substantially leaky scanned. Thus, ATF4's translational control is more complex than originally described, underpinning its key role in diverse biological processes.

Blastocrithidia nonstop mitochondrial genome and its expression are remarkably insulated from nuclear codon reassignment.

The canonical stop codons of the nuclear genome of the trypanosomatid Blastocrithidia nonstop are recoded. Here, we investigated the effect of this recoding on the mitochondrial genome and gene expression. Trypanosomatids possess a single mitochondrion and protein-coding transcripts of this genome require RNA editing in order to generate open reading frames of many transcripts encoded as 'cryptogenes'. Small RNAs that can number in the hundreds direct editing and produce a mitochondrial transcriptome of unusual complexity. We find B. nonstop to have a typical trypanosomatid mitochondrial genetic code, which presumably requires the mitochondrion to disable utilization of the two nucleus-encoded suppressor tRNAs, which appear to be imported into the organelle. Alterations of the protein factors responsible for mRNA editing were also documented, but they have likely originated from sources other than B. nonstop nuclear genome recoding. The population of guide RNAs directing editing is minimal, yet virtually all genes for the plethora of known editing factors are still present. Most intriguingly, despite lacking complex I cryptogene guide RNAs, these cryptogene transcripts are stochastically edited to high levels.

Complete genome sequence of a novel dsRNA virus from the phytopathogenic fungus Fusarium oxysporum.

Fusarium oxysporum is a widespread plant pathogen that causes fusarium wilt and fusarium root rot in many economically significant crops. Here, a novel dsRNA virus tentatively named "Fusarium oxysporum virus 1" (FoV1) was identified in F. oxysporum strain 3S-18. The genome of FoV1 is 2,944 nucleotides (nt) in length and contains two non-overlapping open reading frames (ORF1 and 2). The larger of these, ORF2, encodes an RNA-dependent RNA polymerase (RdRp) of 590 amino acids with a molecular mass of 67.52 kDa. ORF1 encodes a putative nucleocapsid protein consisting of 134 amino acids with a molecular mass of 34.25 kDa. The RdRp domain of FoV1 shares 60.00% to 84.24% sequence identity with non-segmented dsRNA viruses. Phylogenetic analysis further suggested that FoV1 is a new member of the proposed genus "Unirnavirus" accommodating unclassified monopartite dsRNA viruses.

Genome characterization of a novel narnavirus infecting the plant-pathogenic fungus Ustilaginoidea virens.

Mycoviruses are viruses that infect fungi and oomycetes. They are widespread in all major groups of plant-pathogenic fungi and oomycetes. To date, only the full genome of dsRNA mycoviruses and the contigs of positive-sense single-stranded RNA (+ssRNA) mycoviruses have been reported in Ustilaginoidea virens, which is the notorious causal agent of rice false smut (RFS). Here, we report the molecular characterization of a novel +ssRNA mycovirus, Ustilaginoidea virens narnavirus 4 (UvNV4), isolated from U. virens strain Uv418. UvNV4 has a genome of 3,131 nucleotides (nt) and possesses an open reading frame (ORF) predicted to encode an RNA-dependent RNA polymerase (RdRp) of 1,017 amino acids (aa) sequence with a molecular mass of 116.6 kDa. BLASTp analysis revealed that the RdRp showed 50.34% aa sequence identity to that of the previously described Zhangzhou Narna tick virus 1. Phylogenetic analysis indicated that UvNV4 is closely related to members of the family Narnaviridae. Taken together, these results clearly demonstrate that UvNV4 is a novel +ssRNA virus infecting U. virens.

Complete genome sequence of a novel alternavirus isolated from the phytopathogenic fungus Colletotrichum fioriniae.

A novel double-strand RNA (dsRNA) mycovirus, named "Colletotrichum fioriniae alternavirus1" (CfAV1), was isolated from the strain CX7 of Colletotrichum fioriniae, the causal agent of walnut anthracnose. The complete genome of CfAV1 is composed of three dsRNA segments: dsRNA1 (3528 bp), dsRNA2 (2485 bp), and dsRNA3 (2481 bp). The RNA-dependent RNA polymerase (RdRp) is encoded by dsRNA1, while both dsRNA2 and dsRNA3 encode hypothetical proteins. Based on multiple sequence alignments and phylogenetic analysis, CfAV1 is identified as a new member of the family Alternaviridae. This is the first report of an alternavirus that infects the phytopathogenic fungus C. fioriniae.

Complete genome characterization of cacao leafroll virus, a newly described cacao-infecting polerovirus.

The complete genome sequence of cacao leafroll virus (CaLRV; family Solemoviridae, genus Polerovirus) was determined by high-throughput sequencing of total RNA isolated from symptomatic cacao Theobroma cacao L. plants (n = 4). The CaLRV genome sequences ranged from 5,976 to 5,997 nucleotides (nt) in length and contained seven open reading frames (ORFs). Nucleotide and amino acid (aa) sequence comparisons showed that, among selected well-characterized poleroviruses, the CaLRV genome shared the highest nt sequence identity of 62% with that of potato leafroll virus (PLRV, NC_076505). A comparison of the predicted aa sequence of the CaLRV coat protein indicated that cotton leafroll dwarf virus (CLRDV, NC_014545) and melon aphid-borne yellows virus (MABYV, NC_010809) were the closest relatives, sharing 57% aa sequence identity. Bayesian phylogenetic analysis based on complete genome sequences showed that CaLRV grouped with well-characterized poleroviruses that cause diseases of cereal and vegetable crops. During the course of publishing this work, the nearly complete genome sequence of a member of the same polerovirus species, referred to as "cacao polerovirus" (OR605721), with which CaLRV shares 99% nt sequence identity, was reported.

A Precision Therapy Approach for Retinitis Pigmentosa 11 Using Splice-Switching Antisense Oligonucleotides to Restore the Open Reading Frame of PRPF31.

Retinitis pigmentosa 11 is an untreatable, dominantly inherited retinal disease caused by heterozygous mutations in pre-mRNA processing factor 31 PRPF31. The expression level of PRPF31 is linked to incomplete penetrance in affected families; mutation carriers with higher PRPF31 expression can remain asymptomatic. The current study explores an antisense oligonucleotide exon skipping strategy to treat RP11 caused by truncating mutations within PRPF31 exon 12 since it does not appear to encode any domains essential for PRPF31 protein function. Cells derived from a patient carrying a PRPF31 1205C>A nonsense mutation were investigated; PRPF31 transcripts encoded by the 1205C>A allele were undetectable due to nonsense-mediated mRNA decay, resulting in a 46% reduction in PRPF31 mRNA, relative to healthy donor cells. Antisense oligonucleotide-induced skipping of exon 12 rescued the open reading frame with consequent 1.7-fold PRPF31 mRNA upregulation in the RP11 patient fibroblasts. The level of PRPF31 upregulation met the predicted therapeutic threshold of expression inferred in a non-penetrant carrier family member harbouring the same mutation. This study demonstrated increased PRPF31 expression and retention of the nuclear translocation capability for the induced PRPF31 isoform. Future studies should evaluate the function of the induced PRPF31 protein on pre-mRNA splicing in retinal cells to validate the therapeutic approach for amenable RP11-causing mutations.

Molecular and biological characterization of a novel partitivirus from Talaromyces pinophilus.

Talaromyces spp. have a worldwide distribution, are ecologically diverse and have been isolated from numerous different substrates. Talaromyces spp. are considered biotechnologically important due to their ability to produce a range of enzymes and pigments. Talaromyces pinophilus, belonging to genus Talaromyces and family Trichocomaceae, is known for producing several important bioactive metabolites. Here we report the isolation and characterisation of a partitivirus from T. pinophilus which we have nominated Talaromyces pinophilus partitivirus-1 (TpPV-1). TpPV-1 possesses a genome consisting of three double stranded (ds) RNA segments i.e., dsRNAs1-3, 1824 bp, 1638 bp and 1451 bp respectively, which are encapsidated in icosahedral particles 35 nm in diameter. Both dsRNA1 and dsRNA2 contain a single open reading frame (ORF) encoding respectively a 572 amino acid (aa) protein of 65 kDa and a 504 aa protein of 50 kDa. The third segment (dsRNA3) is potentially a satellite RNA. Phylogenetic analysis revealed that the TpPV-1 belongs to the family Partitiviridae in the proposed genus Zetapartitivirus. TpPV-1 infection decreases the mycelial growth rate of the host fungus and alters pigmentation as indicated by time course experiments performed on a range of different solid media comparing virus-infected and virus-free isogenic lines. This is the first report of mycovirus infection in T. pinophilus and may provide insights into understanding the effect of the mycovirus on the production of enzymes and pigments by the host fungus.